www.gusucode.com > 循环自相关函数工具箱源码程序 > matlab代做 修改 程序循环自相关函数工具箱/cyclostationary_toolbox/ex7.m

    clear;
clc;
load('E:\248.mat');  %%读取数据
% speed=1796;    %motor rotational speed(rpm)
sig1=X248_DE_time';
n=1000;
N=512;
t=0:1/12000:(n-1)/12000;
figure(3)
plot(t,sig1(1:n))
x=sig1(1:n).*sig1(1:n);
y=fft(x,N);
Pyy=sqrt(y.*conj(y))/N;
ff=12000*(0:N/2)/N;
figure(4)
plot(ff(1:N/2+1),Pyy(1:N/2+1));